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INTRODUCTION 


This document contains r eference materials, calculations, and trade studies used In the 
analysis and selection of YSTAR components. The document Is organized by each major YSTAR 
system each of which contains material pertinent to that system. Many sections contain simple 
graphs and tables used to make qualitative comparisons of various YSTAR component candidates. 
Equations and/or calculations used for a particular analysis are also included where applicable. 
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Mission Plans 
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aY analysis 


SLR Mission 

This analysis determines the impulse maneuvers required for executing LEO to GEO 
Hohmann transfers for the SLR mission. Two possible methods of transfer are considered for 
this orbit change. The first aV calculation applies to the LEO-GEO transfer involving the 
transfer from LEO to a high-earth orbit followed by a plane change to place the vehicle in 
geosynchronous orbit The total aV required for the method 1 transfer is found to be 
approximately 18,040 ft/sec. The secciw calculation is for a LEO-GEO transfer in which the 
orbital plane change is executed by a two-impulse maneuver within the altitude change. The 
total aV required for the method 2 transfer is found to be approximately 1 4,27 1 ft/sec. 


IYA Mission 

This analysis determines the impulse maneuvers required for executing the 53.27* 
plane change in low-earth orbit ( 28.5* Space Station inclination plus 23.27* Earth inclination 
relative to the eclliptfc plane). The total aV required for the constant orbit plane change is 
approximately 22,504 ft/sec. 

The calculations associated with this analysis assume a low-earth orbit altitude of 230 
miles and a high-earth orbit altitude of 22 ,000 miles. Derivation of equations used in aV 
calculations is given in ( Ref. H ). 
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Table i - Electric Propulsion System Comparisons (**<• 



A&r. high tap 


Mcmaatt 

•IgctrocWMa high lap 
highly davatopad anargy altiaant high thrust 


Occompta siactroda arosion tamp^Bmitad pomaraffidancy highcwrart 

hMrtmd anargy loasaa comdon tSSST* 

Ngh pQMMP "•*****• 

. . .... ... power losses 


Table 2 - Liquid Chemical Rocket Engine Comparisons (*<£ is) 


nocNnoyns v*m » ranpsr vcrsov cn^nvr 

482 sacs Engina langth (ratractad)-1 .02 m 

Propaiant.LH2.LOX (axiandsd)-3.9l m 

Thnjaema*»ae750N Nazzia ax*d ia ma tar- ).98 m 

min-2223 N Arsa ratio- 1300:1 

MR (Off)- 8.0:1 Engina mass-190 kg 


nw era Yvrwwy An^encej um transfer veroae engine: 

bp-486 tecs Engine length (retracted)*1.02 m 

Prape*v*-LH2, LOX (extended)-3.05 m 

' Thrust ma^66750 N Nozzle ex* dteneter* 1.63m 

min-2225 N Area ratio-640:1 

MR(CYF)- 6.0:1 Engine mass-204 kg 


Aerojet Aitanced Ortttd Transfer Vehicle Engine: 

bp ■ 48 2 secs Engine length (retracted)*0.98 m 

f» rops*sr S L H 2.L0X (extended)* 1.96 m 

Thrust max- 13350 N Nozzle ex* ctenetan-0.78 m 

min-445N Ares ratio- 1 200:1 

MR(OF)-6LO:1 Engine meet ■ 57 kg 


Pra* end Whitney RL-10 Category IV: 
bf^470 secs 
PropelanM.H2,LQX 
Thnjet max-66750 N 
min- 1 6668 N 
MR (OF). 6.0:1 


Engine length (retracted)* 1 .45 m 
(extended)*2.90 m 
Nozzle exit demeter*1. 68 m 
Area ratio-401 :1 25 
Engine mess-1 93 kg 
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Table 3 - Dual Fuel / Milled Mode Rocket Engine Characteristics <0 


Pa raw ter 


Mode 1 

(LOX/RP-l/LHj) 


Both 

Modes 


Mode 2 
(lox/lh 2 ) 


Vacuus Thrust (lbf) 20,000. 

Vacuus Delivered 

Specific Inpulse (sec.) 418.6 

Chamber Pressure (psla) 2000. 

Total Propellant Flowrate 

(1b/ sec) 47.8 

Mixture Ratio 4.25 

Nozzle Area Ratio 

Nozzle Exit Olawter (In.) 

Engine Length (in.) 

Total Engine Weight (1b.) 


9915. 

460.6 

1007. 

21.5 

7.0 

400. 

49.3 

6 ! 2/93.2* 

557. 


SINGLE-NODE <yih OTV'S 

OjjHj MIXTURE RATIO -6.1 ^ * 


Oj-COOLED ENGINE; 



DEPLOY P/L - 7926 



SA8A LB 



DRY VT. - 5061 LB 


9478 IB 



» r HIXED-NOOE 01V 

MX. DEPLOY PAYLOAD 
DEPLOY P/L • 9601 LB 

» ■ 



DUAL-FUEL ETJ3. 

NON. DEPLOY PAYLOAD 
7926 LB 


fel-i V 



DRY VT. - 3925 LB 


iynM OTV 

Oj-COOLED ENGINE 
6093 LB 




DRY VT. • 3491 13 




Figure I - Single bs. Mined Mode OTU Comparisons (Ref. b) 
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STRUCTURE 


H 


Thermal Protection 

The selection of the type of radiator was based on survivability, and the membrane area 
per Btu/h output Typically the higher the area required the higher the weight of the device. The 
estimated survivability was based on research. For VSTAR the heat pipes are to be used as thermal 
buses between fluid systems and the choice of these based on the same characteristics as the 
radiator. The data used to determine the ft*2/Btu/h comes from several sources. The Btu/h Is 
converted from Watts and this wattage Is the maximum output of the device under choked flow 
conditions. As can be seen in Table , the best heat pipe for the thermal bus Is the axial ly grooved 
variable-conductance heat pipe, but as a radiator, It's survivability is far to low. As for the 
radiator, the rotating bubble membrane radiator is the superior choice as u is highly survivable 
and does not waste liquid. (Ref. ®i, 


Table Mr. Radiator and Heat Pump Data 


Tvoe 

fTC/Btu/h 

Survivability 

Osmotic Heat Pipe 

0.00293 

medium-low 

Heat Pipe Sandwich Panel 

0.00007 

medium 

Double Walled Artery Heat Pipe 

3.264E-7 

medium -low 

Axially Grooved Variable- 



Conductance Heat Pipe 

1.3I0E-8 

medium 

RBMR 

0.00019 

high 

Liquid Droplet Radiator 

NA* 

medium-high 

Heat Pipe as Radiator (ave) 

0.00075 

medium to low 


* this device sprays the hot liquid directly into space but would weigh the same as the RMBR 
before the liquid was used up. 
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Material Selection 

Table 5 shows the stress allowable and density of each material that was considered. The 
exact cost of each material produced is not known , however , based on research estimated relative 
costs have been determined. Compatibility is known and Is accounted for in the design. The 
material with the smallest densfty/stress is the strongest per unit mass. Metals withstand debris 
and micrometeoroid Impact much better than do composites with epoxy bases, because metals tend 
to buckle before they are punctured and epoxys tend to shatter. 

Therefore a metal Is preferred for debris protection. Careful consideration of this data with the 
micrometeoroid protection results In the choices made for the vehicle components.( Ref. 
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Table 5. Material Properties ( Ref. 2M) 


material 

allowable stress 

density. 

densitv/stress 

cost 


(psl) 

(lbm/1n A 3) 



4 130 Steel 

70000 

0.283 

4.0428E-6 

low 

2024 Aluminum 

42000 

0.100 

2.3809E-6 

very low 

Boron Epoxy 

85333 

0.072 

8.496 IE-7 

medium 

Graphite Epoxy 

100000 

0.057 

5.7000E-7 

medium 

Aluminum Oxide/ 





Aluminum 

42222 

0.100 

2.3684E-6 

low 

Boron Aluminum 





2024 

100000 

0.100 

1.00E-6 

high 

Boron Aluminum 





6061 

96889 

0.100 

1.032E-6 

high 

Graphite Aluminum 





GY70 2024 

37911 

0.100 

2.638E-6 

high 

Graphite Aluminum 





V60054 2024 

43156 

0.100 

2.317E-6 

high 

Graphite Aluminum 





T50 2024 

31111 

0.100 

3.214E-6 

high 



Debris and Micrometeoroid Protection 

BUMPER is a program provided to us by NASA/Marshal to determine the chance of 
penetration of a vessei. BUMPER uses the equations based on some simplifying assumptions and 
your inputs. There are two modes corresponding to single wall and double wall vessels. For double 
wall , inputs are surface area, sheild thickness, sheild stand-off distance, vessel thickness, and 
whether or not MLI is being used. For the single wall surface area, vessel thickness and MU use is 
inputted. The program then returns the percented chance of debris or micrometeoroid protection, 
(Ref. 2.1). 
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Table 1 * EC/LSS Performance Requirements O-** l£ 0 



PARAMETER 

UNITS 

operational ■ 
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Table % - FC/LSS Ruerage Resign Loads •«& 
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Table R - EC/LSS Technology Requirements «<?) 
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FAMILY OF PRIMARY BATTERY CELLS 


MVH9-01DI Hid SUnOH-UVM 



gure ? - Family of Primary Battery Cells 








COST ANALYSIS 
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There ere severe! equations required In cost analysis. These equations are given below for 
each of the subsections, vehicle costs and program management costs. Vehicle costs are In 1 978 dollars 
which must be multlpled by 2 to approximate 1 988 dollars. The management costs are functions of the 
vehicle costs. C***. iO 

Vehicle Costs 

DDT&E 

Structural Cost » (2)( 1.8531 X Structural Weight)®*^ 

Electrical System Cost - (2X0.6000X Electrical System Weight)®* 584 
Propulsion System Cost ■ (2X0.1075X Propulsion System Weight)®* 878 

FH 

Structural Cost » ( 2 )( 0. 4445 )( Structural Weight^- 

Electrical S/stem Cost » (2X0.0422X Electrical System Weight)®* 784 

Propulsion System Cost - (2X0.1 158)( Propulsion System Weight)®* 550 

Management Costs 

Integration, Assembly & Checkout 
DDT&E * 9. 1 208( STE )0* 46 1 
FH - (Total FH)®* 832 

System Test and Evaluation 

DDT&E « 13.1 01 (Total FH) 0 * 672 

FH is zero. It is accounted for in cost of operation. 
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Systems Engineering & Integration 
DDT&E « 0.45 19(DISG) a 876 
FH«0.3750(TotalFH) 0 - 855 

Program Management 
DDT&E - 0.6952(DIS0$) 0,731 
FH«0.3146(FIS) 0 - 798 

Suplementary Equations 

STE » 1 1.232(Totat FH) 0 ’ 672 

DISG * Total DDT&E ♦ Integ., Assem. & Check (DDT&E) + Sys. Test & Eval. (DDT&E) 

DIS6S » STE + Total DDT&E ♦ Integ., Assem. & Check (DDT&E) ♦ Sys. Eng. & Integ. (DDT&E) 
FIS ■ Total FH ♦ Intea . Assem. & Check ( FH) + Sys. Eng. & Eval. ( FH) 
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Optimization Methods 


The following describes the derivation and assumptions made for each of the mathematical 
expressions describing the calculated parameters. The analysis is based on a division of the 
configuration Into the components shown in figure 3 . 

The mass of components in many cases ars tn& -jm of the masses of the subcomponents. 
Their derivations are discussed under the appropriate sections and are not repeated here. 
Calculation of m^^ : The mass of the cabin is constant and not dependent on the 

acceleration loads during flight This is justified by assuming that It will be launched in one piece 
from the space shuttle and will necessarily encounter higher accelerations during launch than 
those experienced for a typical mission. 


m cabin * m crew + m air + m $tructure * m shield1ng + m other 
-600+ 150+ 1140+ 157 + 200 
-2247 

Calculation for mw e support : The mass of the life support system includes the mass of 

the power supply system and is equal to a base mass plus the mass of the consumables ( breathing 
oxygen, water, etc) used per hour by the crew. 

m 1ffe support * f ( * °f c** 8 * » ^end^ 

- 1075+ 1.39 nTgnp 


Calculation of m^go module : The mass of the cargo module is dependent upon the mass of 
structures above it which it will have to accelerate during thrusting periods, as shown in figure 5 
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Figure s 3. Idealized Components for Analysis Purposes 



m cargo module * f ( m cab1n» m 11fe support* m pay* ^nax) 


An approximate function for the mass can be developed by considering figure 
Let A be the cross sectional area of a column of density, r , and maximum allowable stress, 
o, and length 1. The material must withstand the stress generated from mass m and acceleration a. 
The force exerted on the column Is : 

F»ma 


and Is equal : 


F *Ao 

Equating and solving for A gives : 

A»ma/o 

The mass of the column 1$ then : 
m column * A 1 r 

*r1ma/o (1) 

This equation Is based on the assumption that the mass of the column does not contribute 
any to the stresses in It. This is a reasonable assumption if m » matron . 

An extrapolation can be made from the column of Figure 6 to the proposed truss structure 
of the cargo module in Figure s if the additional assumptions are made : 
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a" 



) 


Figure 4 - Stress in o Column 




A 


B 


Figure 5 * Two Possible Tank Configurations 
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1 ) The cross sectional area Is equal to the sum of the cross sectional areas of the 
constituent members of the structure. 

2) The maximum allowable stress In the supporting members Is much lower than the 
compression yield stress to account for lower buckling stresses In a slender column. 

3) The effective density is higher than the actual density to account for the presence of 
cross beams and diagonal members which do not bear the compressive loads. 

Typical r/o values for steel and aluminum are 3.8 x 1 0“6 and 2.9 x 1 0“ 6 [ibm In /lbf] , 
respectively. For Aluminum with the maximum allowable stress to be 1 /1 0 the compressive 
yield stress, the effective density to be 3 times the actual density, and a length of 20 feet, the 
following expression results: 


m support structure * 0.0 1 7 m am®. 


Including an approximate value for the mass of the metiorlte shielding gives : 
m support structure *150 + 0.0 1 7 m 3max 
Where: 


m * m cabin + m hfe support + m pay 


Calculation of m^k ’ Two possible configurations for the arrangement of the fuel tanks 
are shown in figure 5 . 

The twiks in Figures* are of equal radii, but because of different volume requirements 
they are of different lengths. They employ a monocoque structure which resists internal pressure 
and the forces due to acceleration and the mass above it. Figure ?s shows two spherical tanks with 
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a seperate support structue for resisting the acceleration loads. The tanks In Figure are 
efficient in that they provide a dual role In resisting the Internal pressure and the loads due to 
acceleration. However , the nonsherical geometry is Inefficient in providing a light tank for a 
given voluma Additionally, the monocoque skirt between the two tanks will provide an excellent 
path for heat conduction between the two tanks at different temperatures, which if carrying LO 2 
and IH 2 will be over 1 00 degrees ranklne differenca This will complicate thermal control of the 
two tanks and Increase bolloff lossea The arrangement of Figure 78 was chosen because it allows 
the use of spherical tanks which weigh less for a given tank volume, and since they are not actively 
involved in resisting the acceleration loads, they can be more completley thermally isolated. 

The mass of the tank can be expressed In the following form : 

m tank * f ( m fuel» r fuel» m evap» Pfuel* ^mln- r metal) 

The required volume is given by : 

v tank * ( m fuel + m evap ) 1 r fuel 
The radius is given by : 

r tank * ( 3 vtank / 4 pi ) 1/3 

From fluid statics the pressure of a liquid of density, r , under an acceleration, a, at a 
depth of h is : 


p * p 0 + r a h 
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Applying this to a spherical fuel tank.: 


Pbottom * Pvapor + 2 rfuei A^p r^ 

It can be shown that for a constant thrust, the greatest pressure will occur at the bottom 
of the fully filled fuel tank even though the acceleration Is at a minimum. Assuming a thin walled 
pressure vessel the thickness required Is: 

ttank * Pbottom r tank 1 2 o 

The area Is : 

A tank * 4 r tank 2 

The volume of metal In the tank is : 


v metai In tank * A tank ttank 
And the mass of the tank Is: 

m tank * v metal In tank r metal 

Calculation of m^gp: The fuels are cryogenic and are assumed to be continually receiving 
a heat input which holds them at the boiling point for a given pressure. Any teat input is used in 
vaporizing the liquid into its gaseous state, which is immedlatly vented, thus keeping the pressure 
constant. 
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The particular Insulation chosen conducts heat thirty times faster In the direction parallel 
to the layers than In the direction perpendicular to the layers. It will be assumed that the heat 
received from radiation on the Illuminated side of the tank will be conducted to the unlllumlnated 
side allowing an average surface temperature Tout to be used in the below equations. 

It Is assumed that this Is the only source of heat input Into the tanks for the entire 
mission. This is not true however since there will be a radiated and conducted heat transfer during 
the firing of the engines due to the fuel tanks being in close proximity to the hot engine and exhaust 
plume. These calculations may be dependent on many variables, but only the heat due to the 
radiation on the illuminated side of the tank is considered In the following calculation. 

The heat leak rate Into the tank Is given by : 

Odot * k 1nsul A tank (W _t 1n V tinsul 

Then the total heat leak is simply : 

Qtotal * Qdot Tend 

And the mass evaporated is : 

m evap * °total 1 bvapor 

Calculation of m^nsui : The insulation is assumed constant thickness and covers the entire 
spherical tank. 


m 1nsul * r 1n$ul A tank *in$ul 
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Calculation of m^pp^t structure 1 The equation for the mass of the support structure Is 
developed along the same lines as m^rgo module* f t is assumed to be of the form : 

m $upport structure * 200 ♦ 0.0 1 7 m amax 

Where: 


m * m cab1n + m 11fe support + m pay + m tank 
200 - approximate mass of meterolte shielding 

Calculation of m engine • The mass of the engines are assumed to be known for a given 

configuration and are taken from the available literature. 

Calculation of mp rnou t and mtotal : The burnout mass Is simply the sum of all the masses 

calculated above. The total mass Is the sum of the burnout mass and the mass of the fuel. 

Calculation of amax and «mln : The thrust Is assumed constant over the burn time. By 

assuming specific thrust vs. time profiles, the maximum accelerations could be reduced while still 
allowing relatively short burn times, however, here It Is assumed they are constant. 

The accelerations then are simply : 

amax * Thrust / mp rnou t 


amin * Thrust / mtotal 

Calculation of mass ratio : It is assumed their will be no losses due to gravitational and 
drag effects. The mass ratio required then for a given aV and ISP Is : 
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r ■ e aY / g ISP 


Calculation of mfuei : 

From the definition of mass ratio; 

R * ( mbrnout * mfuei ) / mb rnou t 
Solving for mfuei gives: 

mfuei * mbrnout ( R ” 1 ) 

There is one other assumption made In the program that needs to be pointed out. The boiloff occurs 
before the engines Ignite. This results in a higher mass ratio than actually needed since in reality 
the boiloff would occur continually and some of the fuel mass would not be imparted the full 
29,000 ft/sec av. 

Calculation of burn times: 

For a given aY: 

""final * ""initial 8 “ AV/glSP 


The mass of the fuel used is: 

mfuei = ""Initial “ ""final 
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The thrust can be expressed as: 


Thrust * m<jot fuel vexhaust 
Since the values in the above equation are constant- 
Thrust « mfuel g ISP/ 1 
Solving for t 

t » mfuel g ISP / Thrust 
Correctina fnr units: 

t * mfuel g ISP / Thrust ge 
t ■ mfuel ISP / Thrust 

The program listing and program runs are included below. Following this are graphs 
showing data on the liquid hydrogen and liguid oxygen as used in the program. 
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OmvMtftno'* 9«ot,OWA Output 


OXYGEN/HYDROGEN RATIO E . 3 

: 

7 



OXYGEN/RPJ RATIO 

E. 3 

s 

4. 23 



ISP 

C. 3 

X 

420 



ACCELERATION - MAX. 

Eg'a3 

X 

2. 42 



ACCELERATION - MIN. 

Eg'a3 

X 

0. 28 



DELTA V - OBTAINABLE 

Eft/aac3 

X 

29000 



MASS RATIO 

C. 3 

X 

8. 538388 



STUCTUAL COEFF. 

E. 3 

X 

0. 0733 



PAYLOAD COEFFICIENT 

C. 3 

X 

0. 0496 



TIME, ENDURANCE 

C hours 3 

X 

72 



THRUST 

Clbf 3 

X 

30000 



MASSES OF COMPONENTS 

: 





CABIN 

Elba 3 

X 

2247 



LIFE SUPPORT HODUL& 

Elba 3 

X 

1376. 86 



CARGO MODULE 

Elba 3 

X 

504. 8966 



CARGO 

Elba] 

X 

5000 



SUPPORT STRUCTURE 

Elba 3 

X 

642. 4102 



HYDROGEN TANK 

Elba 3 

s 

936. 6639 



OXYGEN TANK 

Elba 3 

: 

448. 6059 



RPJ TANK 

Elba 3 

: 

74. 83677 



ENGINE 

Elba 3 

: 

1000 



TOTAL BURNOUT MASS 

Elba 3 

: 

12392. 82 



MASS OF FUEL 

Elba 3 

: 

93396. 84 



TOTAL MASS 

Elba 3 

: 

105789. 7 






HYDROGEN 

OXYGEN 

RPJ 

TANK - MASS 

Elba3 

t 

936. 66 

448. 61 

74.84 

- MAX. STRESS 

Elbf /in A 23 

: 

35000. 00 

35000. 00 

35000. 00 

- METAL DNSTY 

Elba/ft A 33 

: 

489. 00 

489. 00 

489. 00 

- PRESS. 'BOT. 

Clb£/ln A 23 

: 

20. 14 

21. 67 

21. 28 

- RADIUS 

Eft3 

: 

8.09 

6. 18 

3. 42 

• SURFACE AREA 

Eft A 23 

: 

822. 78 

479. 65 

147. 12 

- THICKNESS 

Ein3 

: 

0. 0279 

0. 0230 

0.0125 

- VOLUME 

Eft A 33 

: 

2219. 22 

987. 80 

167. 79 

- VOLUME. MTL 

Eft A 33 

: 

1.92 

0.92 

0. 15 

PROP. - VAPOR PRS 

Elbf/in A 23 

: 

20. 00 

20.00 

20. 00 

- DENSITY 

Elba/ft A 33 

• 

• 

4.38 

68.64 

95. 11 

- MASS 

Elba 3 

: 

9681.38 

67769. 66 

15945. 80 

- TEMP. 

Erankln»3 

t 

38.38 

166. 41 

58. 33 

- HEAT. VAPRZ 

Ebtu/lbm] 

X 

188. 16 

91.62 

100. 00 

SURFACE TEMPERATUE 

Crankln*] 

• 

• 

600. 00 

600. 00 

600. 00 

INSULATION - RHO Elba/ft A 2-in3 

: 

0.20 

0.20 

0. 20 

- MASS 

Elba 3 

: 

41. 14 

23.98 

7. 36 

- THCKNS 

Eln3 

• 

• 

0.25 

0.25 

0. 25 

-K Ebtu-in/hr-ft A 2J 

• 

• 

0. 000056 

0. 000056 

0. 000056 

HEAT LEAK - RATE 

Ebtu/hour 3 

: 

103. 51 

46. 59 

17. 85 

- TOT 

Cbtu3 

• 

• 

7452. 54 

3354. 22 

1285. 23 

MASS EVAPORATED 

Elba 3 

: 

39.61 

36.61 

12. 85 


36 



OXYGEN/HYDROGEN RATIO C. 3 


ISP 

C. 3 

ACCELERATION - MAX. 

Cg's3 

ACCELERATION - WIN. 

Cg's3 

DELTA V - OBTAINABLE 

Cf t/sac3 

HASS RATIO 

C. 3 

STUCTUAL COEFF. 

C. ] 

PAYLOAD COEFFICIENT 

C. 3 

TIHE, ENDURANCE 

C hours 3 

THRUST 

(lbf 3 

HASSES OF COHPONENTS 

• 

• 

CABIN 

Clbm3 

LIFE SUPPORT HODULE 

Clbm 3 

CARGO HODULE 

Clbm 3 

CARGO 

Clbm3 

SUPPORT STRUCTURE 

Clbm 3 

HYDROGEN TANK 

Clbm3 

OXYGEN TANK 

Clbm 3 

ENGINE 

Clbm 3 

TOTAL BURNOUT HASS 

Clbm3 

HASS OF FUEL 

Clbm3 

TOTAL HASS 

Clbm] 


TANK 

- 

HASS 

Clbm 3 


• 

HAX. STRESS 

Clbf /ln A 23 


• 

HETAL DNSTY 

Clbm/ft A 33 


- 

PRESS. , BOT. 

Clbf/in A 23 


- 

RADIUS 

Cft 3 


- 

SURFACE AREA 

Cft A 23 


- 

THICKNESS 

[in 3 


- 

VOLUHE 

Cft A 33 


- 

VOLUME, MTL 

If t A 33 

PROP. 

- 

VAPOR PRS 

Clbf /in A 23 


- 

DENSITY 

Clbm/ft A 33 


- 

HASS 

Clbm 3 


- 

TEHP. 

Crankins3 


- 

HEAT. VAPRZ 

Cbtu/lbm3 

SURFACE 

TEMPERATUE 

Crankin*3 

INSULATION - RHO Clbm/ft A 2-in 3 



- HASS 

Clbm 3 



- THCKNS 

Cin3 



-K Cbtu-in/hr-ft A 23 

HEAT 

LEAK - RATE 

Cbtu/hour 3 



- TOT 

Cbtu3 

HASS 

EVAPORATED 

Clbm 3 


: 7 
: 460 
t 2.45 
: 0.35 
: 29000 
: 7.084234 
2 0.0887 
: 0.0611 
: 72 
: 30000 


2 2247 
: 1376.86 
2 508.9224 
2 5000 
2 641.6741 
2 902.9588 
t 438.6012 
: 1000 
: 12253.81 
: 74555.06 
2 86808.88 


HYDROGEN 

OXYGEN 

902. 96 

438. 60 

35000. 00 

35000. 00 

489. 00 

489. 00 

20. 17 

22.01 

7.99 

6. 10 

802. 17 

467.62 

0. 0276 

0. 0230 

2136. 36 

950. 87 

1.85 

0. 90 

20.00 

20. OO 

4.38 

68. 64 

9319. 38 

65235. 68 

38.38 

166. 41 

188. 16 

91. 62 

600. 00 

600. 00 

0. 20 

0. 20 

40. 11 

23. 38 

0.25 

0. 25 

0. 000056 

0. 000056 

100.91 

45. 42 

7265. 84 

3270. 10 

38.62 

35. 69 
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OXYGEN/HYDROGEN RATIO C. 3 


ISP 

C. ] 

ACCELERATION - WAX. 

Cg'a] 

ACCELERATION - WIN. 

Cg'a] 

DELTA V - OBTAINABLE 

Cft/sec] 

HASS RATIO 

C. 1 

STUCTUAL COEFF. 

C. 3 

PAYLOAD COEFFICIENT 

C. 3 

TIME, ENDURANCE 

C hours 3 

THRUST 

Clbf 3 

MASSES OF COMPONENTS 

: 

CABIN 

Clbm3 

LIFE SUPPORT MODULE 

Clbra] 

CARGO MODULE 

Clbra3 

CARGO 

Clbra] 

SUPPORT STRUCTURE 

Clbra] 

HYDROGEN TANK 

Clbra] 

OXYGEN TANK 

Clbra] 

RPJ TANK 

Clbra] 

ENGINE 

Clbra] 

TOTAL BURNOUT MASS 

Clbra] 

MASS OF FUEL 

Clbra] 

TOTAL MASS 

Clbra] 


TANK 

• 

MASS 

Clbra] 


mm 

MAX. STRESS 

Clb£/in A 2] 


- 

METAL DNSTY 

C lbm/ft A 31 


- 

PRESS. , BOT. 

Clb£/in A 2] 


- 

RADIUS 

tft] 


- 

SURFACE AREA 

Cft A 23 


- 

THICKNESS 

tin] 


- 

VOLUME 

C ft A 3 ] 


- 

VOLUME, MTL 

Cft A 33 

PROP. 

- 

VAPOR PRS 

C lb£/in A 23 


- 

DENSITY 

C lbrn/ft A 33 


- 

MASS 

Clbra] 


- 

TEMP. 

Crankine] 


- 

HEAT. VAPRZ 

Cbtu/lbm] 

SURFACE 

TEMPERATUE 

Crankine] 

INSULATION - RHO C Ibm/ft A 2-in 1 



- MASS 

Clbra] 



- THCKNS 

Cin] 



-K Cbtu-in/hr-ft A 23 

HEAT 

LEAK - RATE 

Cbtu/hour ] 



- TOT 

C btu ] 

MASS 

EVAPORATED 

Clbra] 


6 

470 
2.44 
0. 36 

28999. 99 
6. 795188 
0. 0930 
0. 0636 
72 

• 30000 


j 2247 


: 1376.86 


: 507. 1932 


: 5000 


: 641.9901 


: 988. 1334 


: 411.8053 


: 3. 70943E-08 


: 1000 


: 12313. 14 


s 71356.95 


s 83670.09 


HYDROGEN 

OXYGEN 

988. 13 

411. 81 

35000. 00 

35000. 00 

489. 00 

489. 00 

20. 18 

22. 04 

8.23 

5.97 

851. 53 

447. 96 

0. 0285 

0. 0226 

2336. 53 

891. 52 

2. 02 

O. 84 

20.00 

20. OO 

4. 38 

68. 64 

10193.85 

61163. 10 

38. 38 

166. 41 

188. 16 

91. 62 

600. 00 

600. 00 

0. 20 

0. 20 

42. 58 

22. 40 

0. 25 

0. 25 

0. 000056 

0. 000056 

107. 12 

43. 51 

7712. 90 

3132. 57 

40. 99 

34. 19 


38 



OXYGEN/HYDROGEN RATIO E. 3 

t 

6 


ISP 

C. 3 

2 

486 


ACCELERATION - MAX. 

Eg'aJ 

X 

2. 46 


ACCELERATION - MIN. 

Eg'aJ 

X 

0.39 


DELTA Y - OBTAINABLE Eft/aac3 

X 

29000 


HASS RATIO 

C. 3 

: 

6. 379734 


STUCTUAL COEFF. 

C. 3 

: 

0. 0988 


PAYLOAD COEFFICIENT 

C. 3 

X 

0. 0687 


TIHE, ENDURANCE 

[hours 3 

X 

72 


THRUST 

Clbf 3 

X 

30000 


HASSES OF COHPONENTS 

t 




CABIN 

Elbal 

X 

2247 


LIFE SUPPORT HODULE 

ElbaJ 

X 

1376. 86 


CARGO HODULE 

ElbaJ 

X 

310. 3693 


CARGO 

Clbm3 

X 

3000 


SUPPORT STRUCTURE 

Elba! 

X 

641. 3731 


HYDROGEN TANK 

ElbaJ 

X 

909. 1804 


OXYGEN TANK 

Elba 3 

X 

380. 309 


ENGINE 

Elba 3 

X 

1000 


TOTAL BURNOUT HASS 

Elba 3 

X 

12197. 84 


HASS OF FUEL 

Elba 3 

X 

65621.38 


TOTAL HASS 

Elba 3 

X 

77819. 22 





HYDROGEN 

OXYGEN 

TANK - HASS 

Elba 3 t 


909. 18 

380. 31 

- HAX. STRESS 

Elbf/in A 23 * 


35000. 00 

33000. 00 

- HETAL DNSTY 

Elba/ft A 33 : 


489. 00 

489. 00 

- PRESS. , BOT. 

Elbf/in A 2J: 


20. 19 

22. 13 

• RADIUS 

Eft3 : 


8.01 

3.81 

- SURFACE AREA 

Eft A 2J : 


805. 32 

423. 62 

- THICKNESS 

Eln 3 t 


0. 0277 

0. 0220 

- VOLUHE 

Eft A 33 t 


2148. 97 

819. 87 

- VOLUHE, HTL 

Eft A 33 : 


1.86 

0.78 

PROP. - VAPOR PRS 

Elbf/in A 23: 


20.00 

20.00 

- DENSITY 

Elba/ft A 33 : 


4.38 

68. 64 

- HASS 

Elba 3 : 


9374. 48 

56246. 90 

- TEHP. 

Erankina J : 


38. 38 

166. 41 

- HEAT. VAPRZ 

Ebtu/lbaJ : 


188. 16 

91.62 

SURFACE TEHPERATUE 

Erankina 3 : 


600. 00 

600. 00 

INSULATION - RHO Elba/ft A 2-inJ s 


0. 20 

0. 20 

- HASS 

Elba 3 : 


40. 27 

21. 18 

- THCKNS 

Eln 3 : 


0.23 

0. 25 

-K Ebtu-in/hr-ft A 2J : 


0. 000036 

0. 000056 

HEAT LEAK - RATE 

Ebtu/hour 3 : 


101.31 

41. 14 

- TOT 

Ebtu3 : 


7294. 41 

2962. 40 

HASS EVAPORATED 

Elba 3 : 


38.77 

32.33 
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V 5TAR. O 


10 00SUB 1000 'INITIALIZE VARIABLES 
20 CLS*INPUT a OPTION 1 OR 2 •|IJOB ... 

30 ON I JOB GOTO 40,210 
. 40 OPEN a O a ', #1, "R0CXET1. DAT* 

30 FOR CONFIG-1 TO 6 * 

60 READ OXH, OXRP, ISP 
70 TEND-72 
« 80 THRUST -30000 

90 TINSULd > 23 

100 TINSUL(2)-.2S 
110 TINSUL(3) •» 23 
120 PVAPOR < 1 > -20 
130 P VAPOR ( 2 > -20 
140 PVAPOR ( 3 ) -20 
130 GOSUB 2000 
160 GOSUB 3000 'GET MFUEL 

170 PRINTfl, OXH, OXRP, ISP, TEND, THRUST, TINSUL ( 1 ) , TINSUL ( 2 > , TINSUL ( 3 ) , PVAPOR ( 1 ) , 
POR < 2 ) , PVAPOR ( 3 ) , MFUEL 

173 PRINT OXH, OXRP, ISP, TEND, THRUST, TINSUL ( 1 ) , TINSUL < 2 > , TINSUL ( 3 ) , PVAPOR ( 1 ) , PV 

R ( 2 ) , PVAPOR ( 3 > ►'FUEL 

180 NEXT CONFIG 

190 CLOSE# 1 

200 STOP 

210 GOSUB 7000 ' SORT DATA 
220 OPEN*O a , #1, a R0CKET3. DAT* 

230 OPEN a I a , #2, a R0CKET2. DAT a 

240 INPUT#2, OXH, OXRP, ISP, TEND, THRUST, TINSUL < 1 ) , TINSUL ( 2 ) , TINSUL ( 3 > , PVAPOR ( 1 ) 

POR < 2 > , PVAPOR ( 3 ) , HFUEL 

230 GOSUB 2000 ' GET FUEL PROPERTIES 

260 GOSUB 3000 ' GET PARAMETERS 

270 GOSUB 4000 * PRINT PARAMETERS 

280 IF EOF ( 2) -0 THEN 240 

290 CLOSE 

300 STOP 


/ 
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1000 ' ••••••••••#••••••••• INITIALIZE CONSTANTS 

iS5t If!' RVAP0R <3) ' MEVAP <3) ' NAXSTRS < 3 > , MTANK < 3 ) , PA < 3 > , PTEHP < 3 ) , P VAPOR 
^QDOTO), QT0T(3), RADIUS < 3 ) , RFUEL ( 3 ) , RINSUL(3>, RHETAL(3), TANK*(3), TINSULO), TOl 

101G DELTA V0-290001 
1020 KINSULU)-. 000036 
1020 KINSUL<2) ■. 000036 
1040 KliiSULO)-. 000036 


10S0 HAXSTRS(l) -330001 
1090 NAXSTRS (27 -330001 
1100 NAXSTRS ( 3 >-330001 
1110 HPAY-3000 
1120 PI-3.1413927* 

1130 RHETAL ( 1 ) -489 
1140 RHETAL (2) -469 
1130 RHETAL ( 3 ) *469 
1160 R1NSUL1 1 ) -. 2 
1170 RINSULC2) -. 2 
1160 RXNSUL(3) 2 
1190 TANKA(l) ••HYDROGEN* 

1200 TANK*< 2) ••OXYGEN" 

1210 TANK813) -*RPJ* 

1220 TOUT ( 1 ) -600 
1230 TOUT (2) -600 
1240 TOUT ( 3 ) -600 
1230 DIN N(3, 3) 

1260 FOR X-l TO 3sF0R Y«1 TO 3 : READ N ( X, Y > * NEXT Y : NEXT X 
1270 DATA 8.10.6 
1280 DATA 10.3,2 
1290 DATA 3,3,2 

1300 DIN THV APOR ( 3, 6 ) , HVAP0R1 (3,6), PFUEL (3,10), TPFUEL (3,10), RFUEL1 (3,10), TRFUE 
3,10) 

1310 FOR X-l TO 3 

1320 FOR Y-l TO N(X, 1):READ PFUEL(X, Y>, TPFUEL (X, Y) t NEXT Y 
1330 FOR Y-l TO N(X,2):READ RFUEL1 ( X, Y ) , TRFUEL ( X, Y ) t NEXT Y 
1340 FOR Y-l TO N ( X, 3 ) t READ THVAPOR ( X, Y ) , HVAPQR1 ( X, Y ) : NEXT Y 
1330 NEXT X 
1360 RETURN 

1370 DATA 1. 044, 23. 07, 1. 102, 23. 27, 2. 939, 28. 67, 6. 689, 32. 27, 13. 066, 33. 67, 22. 992, 
.47,29.446,41.27,33.903,42.37 

1380 DATA 4. 71, 26. 87.,4. 61, 30. 47, 4. 46, 34. 07, 4. 43, 37. 67, 4. 16, 41. 27, 3. 99, 43. 27, 3. 
, 48. 47, 3. 3, 32. 07, 3. 14, S3. 67, 2. 22. 39. 27 

1390 DATA 32. 27, 898. 83, 33. 87, 860. 16, 39. 79, 804. 12, 43. 4, 744. 19, 47, 674. 61, 30. 6, 5« 


22 

1400 DATA . 027, 98. 67, . 033, 116. 67, 2. 102, 134. 67, 8. 223, 132. 67, 23. 633, 170. 67, 76. 7 J 

197. 67, 196. 692, 224. 67, 322. 342, 242. 67, SOI. 723, 260. 87, 611. 814, 269. 67 

1410 DATA 83. 63, 97. 77, 71. 2, 162. 77, . 2766, 263. 67 

1420 DATA 10,91.62,300,91.62 

1430 DATA 10,33,11,33,12,33 

1440 DATA 73,38,73,42,73,46 

1430 DATA 20, 100, 70, 100 

1460 DATA 7,2.21,460 

1470 DATA 7,1000000000000000000,460 

1480 DATA 7, 4. 23, 420 ' 

1490 DATA 6,1000000000000000000,492 
1300 DATA 6, 1000000000000000000, 486 
1310 DATA 6, 1000000000000000000, 470 
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2000 »*****•**#**.•**•*•* GET FUEL DENSITIES 

2010 FOR X»1 TO 3 •'.•••- 

2020 FOR Y»1 TO N(X , It / . 

2030 IF PVAPORtXXPFUEUX, Y> THEN 2030 
2040 NEXT Y 

2030 PTEMP ( X > -tPFUEL < X, Y-l > ♦ < TPFUEL < X, Y ) -TPFUEL < X, Y-l ) ) * < PVAPOR < X ) -PFUEL < X, Y -1 ) 
/ < PFUEL < X, Y ) -PFUEL (X, Y-l) > 

2060 FOR Y«1 TO N<X, 2) 

2070 IF PTEMP ( X X TRFUEL <X,Y> THEN 2090 
2060 NEXT Y ' 

2090 RFUEL ( X > •RFUEL1 ( X, Y-l > * < RFUEL1 < X, Y ) -RFUEL1 < X, Y-l > > * < PTEMP < X ) -TRFUEL <X, Y-l > 
/ ( TRFUEL ( X, Y ) -TRFUEL ( X, Y-l ) ) 

2100 FOR Y»1 TO N<X,3) 

2110 IF PTEMP (XX THVAPOR<X,Y> THEN 2130 
2120 NEXT Y 

2130 H VAPOR < X > -HVAPORl ( X, Y-l ) ♦ < HVAP0R1 ( X, Y ) -HVAP0R1 ( X, Y-l ) ) • < PTEMP < X ) -THVAPOR < X 
Y-l > > / < THVAPOR ( X, Y ) -THVAPOR ( X, Y-l ) ) 

2140 NEXT X 

2130 H VAPOR ( 1 ) -HVAPOR ( 1 ) /RFUEL ( 1 ) 

2160 RETURN 
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3000 '•*♦**•**♦•****••*•♦ FIND CABIN HASS 
' 3010 HBRNOUT-O 

3020 HCA*BIN-600-130-137-1140-100-100 
3030 HBRNOUT-HCABIN 

3040 '•***•♦•*••*•♦••***• FIND LIFE SUPPORT SYSTEM HASS 

3030 HLIFSUP-1073* (33. 34/24) *3*TEND 

3060 HBRNOUT-HBRNOUT-MLIFSUP-MPAY 

3070 '**•**•*•*****••***• FIND CARGO MODULE HASS 

3080 MCRGH0D"130*.017*HBRN0UT*AMAX 

3090 MBRNOUT* MBRNQUT-MCRGMOD 

3100 '**••*••**•****♦*♦••* FIND TANK MASSES 

3110 HFUEL ( 1 ) "MFUEL/ ( 1-OXH-OXH/OXRP ) 

3120 HFUEL < 2 ) * HFUEL ( 1 ) *OXH 
3130 HFUEL ( 3 ) "HFUEL ( 2 ) /OXRP 
3140 FOR X"1 TO 3 
3130 MEVAP»0 

3160 VTANK < X > • < HFUEL ( X ) -MEVAP ( X ) ) /RFUEL < X ) 

3170 RADIUS < X ) • < 3* VTANK < X > /PI/4 > A < 1/3 > 

3180 AREA(X) *4*PI*RADIUS(X) *2 

3190 PA < X ) "PVAPOR < X ) * < RFUEL < X ) / 12“ 3 ) • AMIN* < 2*RADIUS < X ) • 12 ) 
3200 TA ( X ) "PA ( X ) * RADIUS ( X ) • 12/2/ HAXSTRS ( X ) 

3210 T<X)«TA<X) 

3220 VMETAL ( X ) "AREA <X)*T(X)/12 
3230 HTANK ( X ) * VMETAL ( X ) • RHETAL ( X ) 

3240 ODOT(X) "KINSUL(X) *AREA (X) * <TQUT<X> -PTEMP<X> ) /TINSULC X ) 
3230 QTOT ( X ) "QDOT ( X ) "TEND 

3260 'PRINT USING"####. ####-jMEVAP, MEVAP(X) 

3280 MEVAP < X ) "QTOT ( X ) /HV APQR l X ) 

3290 IF ABS(MEVAP-HEVAP<X) )>.01 THEN MEVAP-MEVAP<X) sGOTO 3160 
3300 HINSUL ( X ) "RINSUL ( X ) * AREA < X ) »TINSUL ( X > 

3310 MBRNOUT-MBRNOUT-MTANK ( X ) -HINSUL ( X ) -HEVAP ( X ) 

3320 NEXT X 

3330 '•**.**••*♦*.•••***♦• FIND SUPPORT STRUCTURE MASS 

3340 MSUPSTR • 200 - . 017*HBRN0UT»AMAX 

3330 H BRNOUT > M BRNQUT * MSUPSTR 

3360 '•♦*♦•****•***♦**•*•* FIND ENGINE MASS 

3370 HENGINE" 1000 

3380 '*•••*♦*••***♦••**•*• FIND BURNOUT HASS 
3390 MBRNOUT*MBRNOUT-MENGINE 
3400 HTOTAL * MBRNOUT * MFUEL 
3410 AMAX-THRUST/MBRNOUT 
3420 AMIN-THRUST/MTOTAL 

3430 'PRINT USING”#####. ####* ;AMAX, AMIN. AC 
3440 IF ABS ( AMAX-AC > >. 01 THEN AC*AMAX>GOTO 3000 
3430 MRATI0*2. 7182818# A ( DELTAVO/32. 2/ISP ) 

3460 MFUEL1" MBRNOUT* ( MRATIO-1 ) 

3470 'PRINT USING”##########. ##• ;MFUEL1 ; MFUEL ;MRATIO 
3480 IF ABS< HFUEL1 -HFUEL )>. 01 THEN HFUEL" MFUEL1 sGOTO 3000 
3490 HRATIQ-MTOTAL/ MBRNOUT 
3300 DELTAV-32. 2*ISP*L0G( MRATIO) 

3310 EPSILON" ( HBRNOUT-MPAY ) / ( HTOTAL-MPAY ) 

3320 LAMBDA«HPAY/< HTOTAL-MPAY) 

3330 RETURN 
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4120 

PRXNT#1, 

• CABIN 

ClbsJ 

< • ; MCABIN 

4130 

PRXNT#1, 
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t " ; MLXFSUP 

4140 

PRXNT#1, 

" CARGO MODULE 

ClbsJ 

j ■ ; HCRGHOD 

4190 

PRINT#1, 

• CARGO 

ClbsJ 

s " >HPAY 

4160 

PRINT#1, 

” SUPPORT STRUCTURE 

ClbsJ * 

* "jMSUPSTR 

4170 

PRINT#1, 

” HYDROGEN TANK 

ClbsJ 

J " ; MTANK ( 1 ) 

4180 

PRINT#1, 

• OXYGEN TANK 

ClbsJ 

s " ; MTANK < 2 > 

4190 

PRINT#!, 

" RPJ TANK 

ClbsJ 

1 " {MTANK (3* 

4200 

PRINT#!, 

" ENGINE 

ClbsJ 

: " j MENGINE 

4210 

PRINT#!, 

"TOTAL BURNOUT HASS 

ClbsJ 

: " ; HBRNOUT 

4220 

PRINT#1, 

"HASS OF FUEL 

ClbsJ 

j • ; HFUEL 

4230 

PRINT*!, 

"TOTAL MASS 

ClbsJ 

: • ; HTOTAL 
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, • , i PRINT# 1# USING •##########- #*" i 
i *> »PRINT#1, USING" #######***• #*" »■ 
, • , *PRINT#1, USING" ######*#• ***** 
t • f tPRINT#!, USING"#####****'- ** * J 
, • i jpftlNT#!, USING"##********- **" J 


SCKW PRIHT#i,TAB(2S),"TANK # ’’ TA ***^ 2 SING -########*#.*♦' 

9010 PRINT#1, "TANK - HASS !; AB ’ J / 

ANK< 1), HTANK<2)» STANK < 3 K £ Ib*/±n*21 i • i *PRIHT#1, USING*#######**'- » 

3020 PRINT# 1, " “ MAX " STRESS uw '* l ‘ 

SS‘re'xS«!” S<2> ' “ X «^ 3 B«srr cu»/«-ai * * > •«««"• 

SSo^SS#^^’' . BOT. tibi/in-21 . • » .PRINT#!. USIHS •».....»•** > 

. Radius ««i 

IUS11), RADIUS <2 ), RADIUS (3) £ft*2J 

SOSO PRINT#1, * “ SURFACE AREA £« 2J 

A ( 1 ) • AREA < 2 > # AREA ( 3 ) 

5070 PRINT#!# • “ THICKNESS iinl 

<1># T(2)*T(3> volume 

5080 PRINT#i# • • VOLUHE ** 

NK<1>,VTANK12),VTANK13J 

3090 PRINT#!, ■ - VOLUNE, HTL Ctt 

St^PRwS^ PRB . • » .PRIHT.l. USING-........**- **• “ 

* 2> ' PV *-“^SIT» Clb./«*3! .PRIHT.l, USING -......*•••- **' 

EL(1),RFUEL(2),RFUEL(3) , • , »PRINT#1, USING "####*##***- ** * ; 1 

3120 PRINT#!, " - «Aas lAO " J 

«3i‘p"lN«if " FUEL - 3 iENP. irukim j > • I .PRINT#!. USING-....#.#*##. ** ' > f 

Mio’pSmCif •' "“‘‘-^EAT. V.PR2 CM»/lb.l • -PRIHT.l. USING-.........#- 

S130 1 pRINT»l^ -SURPACE°TNNPtRA TUE Cr.nNin.J * * » .PRINT.l. USING -.....*.»*»• ## " » ^ 

;iii* ss;;i:'iS2i2iviu« - r-g cib./«-2-i»i . • » ^^.1. usihu....... — - 

SUH 1 /, RINSUV.12J. RINSULiJ) ' . st »RINV#l, USING" #*»«*#*##*- *#'* ft 

3170 PRINT#1, " _* LAO 

suni),ni**sun2>.«iN»unji * • 1 « print#!, using •####***##*• ** " * ^ 

S1S0 PRINT#1, " " THuhHi> ixnj 

SU ^ A ^Ti^ l 'i 2,,TINSUl ' V3 -K CStu-xn/hr-it. A 23 * " J i PRINT#!, USING"######* ###***" J ^ 
5 1 SO « R in 1 #•( 

So^PrS ?^ 1 ^^"^ 41 rate Cbtu/ Hour 1 . - ; . PRINT.!, USING — — - 

T 1 1 > , ODOT < 2 * , QDOT < 3 ) t • ; s PRINT#1, USING* ####*##***• **";•• 

3210 PRINT#!, " * TOT lDlu -’ 

Til), QTOT 12), OTOT ( 3 ) c ibB j PRINT#! , USING " *#*#*####*. ** " ) « 

3220 PRINX*1, " BA!s! » EVAPURATED llo " J 
APll), HEVAP12), NEVAP13) 
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6000 HANV6 < 1 ) - * INJECT HOHHAN TRANSFER/* 

6010 HANV6< 2) -"CIRCULARIZE '*. 

6020 HANV«<3>«*D£-CIRCULARIZE >• . 

6030 HANV«<4> -'CIRCULARIZE • 

6040 OELTAV ( 1 ) -6329. 87 
6030 0ELTAV(2>-3741. 13 
6060 DELTAV(3)-3741. 13 
• 6070 OELTAV < 4 > *6329. 67 
6060 T-HFUEL-ISP/THRUST/60 
6090 PRINT#1, "TIME OF BURN CminutM] r*;T 
6100 H1N1T-HT0TAL 
6110 PRINT* FOR X-l TO 4 

6120 HFINAL-M INIT-2. 72* ( -DELTA V ( X ) /32. 2/ ISP > 
6130 HFUEL-HINIT-HFINAL 
6140 T-HFUEL-ISP/THRUST/60 

6130 PRINT#!,* - * j HAN V* (X) j*i*jT, DELTA V ( X > 

6160 HINIT-HINIT-HFUEL 
6170 NEXT X 

6160 FOR X-l TO 13:PR1NT#1, *”:N£XT X 
6190 RETURN 
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7010 OIK J 4 ; gJjcKETl.OAT- . 

&£&****'■ , 
WWWfP* - 

7130 A<Y>'-Al*®« Y,- f, W 



7150 PRIHT#2, A 

V$> Tr*n<-* ™ EH 70*0 
7180 CLOSE’ 

7190 RETURN 
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6000 OPEN'O*, #1, -FUEL. DAT* / 

6010 FOR X«1 TO 3 ' 

6020 FOR Y*i TO N(X # 1) »PRINT#1, FFUELOC, Y),TPFUEL(X, Y>.iNEXT Y 
6030 FOR Y«1 TO N < X, 2 ) t PRINT#1, RFUEL1 ( X, Y > , TRFUEL < X» Y > I NEXT Y 
6040 FOR M TO N<X, 3) *PRXNT#1, THVAPOR<X, Y>, HVAP0R1 ‘.X, Y) jHEXT Y 
6090 NEXT X 
6060 CL0SE#1 
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Temperature at boiling point vs. Pressure 
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Pressure (PSI) 




OXYGEN 

Density of liquid at boiling point vs. Temperature 
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